Quantum and Classical Dynamics of a BEC in a Large-Period Optical Lattice 
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We experimentally investigate diflraction ol a *^Rb Bose-Einstein condensate from a ID optical 
lattice. We use a range of lattice periods and timescales, including those beyond the Raman-Nath 
limit. We compare the results to quantum mechanical and classical simulations, with quantitative 
and qualitative agreement, respectively. The classical simulation predicts that the envelope of the 
time-evolving diffraction pattern is shaped by caustics: singularities in the phase space density of 
classical trajectories. This behavior becomes increasingly clear as the lattice period grows. 
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Introduction 

Modern atom optics [TJ [5] , in particular the diffraction 
of atoms by standing waves of light H] provides a dra- 
matic demonstration of the wave nature of atoms. The 
advent of Bose-Einstein condensates (BEC), with their 
extremely narrow momentum distribution, has made re- 
solving diffraction components straightforward in these 
systems [S]. Most atom diffraction experiments with 
BECs focused on the regime where the diffracting stand- 
ing wave has relatively few "bound" states [6] or where 
the diffraction produces relatively few diffraction orders. 
Under such conditions, the wave nature of the atoms is es- 
sential for describing the behavior of the system. By con- 
trast, when the optical potential has many bound states 
the quantum system can be well described by classical 
particle trajectories. This latter regime has been theo- 
retically investigated with such trajectories [7] , as well as 
from a matter wave perspective [5J [3] . 

In 1994, Janicke and Wilkens calculated the diffraction 
of cold atoms from a lattice potential and predicted a 
dramatic collapse and revival of the number of diffracted 
orders as the atoms coherently evolve in the lattice [S]. 
Here we present the first experimental observation of this 
collapse and revival in both the quantum and classical 
regimes. The envelope of the time-evolving diffraction 
pattern for a sinusoidal potential is dominated by caus- 
tics: singularities in the phase space density of classical 
trajectories |7]. We investigate the quantum and classical 
regimes of atom wave diffraction by applying an optical 
standing wave (an optical lattice) to a BEC, and mea- 
suring the time evolution of the momentum distribution 
for a range of lattice periods. 

While the phenomenon of atom diffraction is quan- 
tum mechanical, some early experiments used a classi- 
cal trajectory approach to describe the observed channel- 
ing [Tir and focussing [TTl [T^] of atoms by optical stand- 
ing waves. These experiments used thermal beams of 
atoms with a large momentum spread and measured the 
atomic position distribution within the standing wave. 



Another early experiment used an optical standing wave 
to diffract laser-cooled atoms and observed the growth 
and subsequent collapse of the width of the diffraction 
pattern in the quantum regime |13| . Later, the collapse 
and revival of a few diffraction orders was observed in 
the diffraction of a BEC from a lattice with only a few 
bound states |14j . Here we extend this earlier work by 
measuring the time-evolution of a BEC's momentum dis- 
tribution in lattice potentials for a range of lattice periods 
(see Fig. 1). We observe several collapses and revivals in 
both quantum and classical regimes. We compare our re- 
sults to the predictions of a single-particle quantum sim- 
ulation and find excellent agreement. We also develop a 
classical model which reproduces essential features of our 
data and provides physical insight into the evolution of 
the momentum distribution. 

Experimental Procedure 

We create a static ID optical lattice at the intersection 
of two laser beams. In the plane- wave approximation, the 
electric field of each beam is E{f, t) = eEoe^^'^'^'^^*''' + c.c. 
where e is the polarization vector. The combined electric 
field for both beams creates an optical potential for the 
atoms given by 

U{z) ^Uasin^i^Lz). (1) 

For a two-level system Uq = fiVl^/d, and kl = \ki — 
^21/2 — n/d is one-half of the magnitude of the recipro- 
cal lattice vector. A/ is the atomic mass, d is the lat- 
tice period, f2o is the resonant, single-beam Rabi fre- 
quency, and 5 < is the detuning of the laser from 
atomic resonance (In this experiment, C/q is nominally 
30 Eji where E-^i = h^k^/2M is the single-photon recoil 
energy.). We vary kl by changing the angle between ki 
and k2 and define z parallel to fci — ^2. Near the energy 
minimum of each lattice site, the potential is nearly har- 
monic, Ui,o{z) « where = 2Uan^/Md'^. 
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We choose the detuning |i5| ^ (flo, T) (where F is the nat- 
ural linewidth of the atomic transition) so spontaneous 
emission is neghgible for our experiment durations. 

We produce a nearly pure BEC with iVtot = 5 x 10'' 
to 14 X 10^ atoms in the {F,mp) = (1,-1) hyperfine 
state of ^^Rb [15 . We use a loffc-Pritchard trap with an 
oscillation frequency of — 8.2 Hz in the weak direction 
and v^ — Uy — 24 Hz in the tight directions. The optical 
lattice is suddenly applied to the magnetically-trapped 
atoms for a time T^vlIsc- The lattice periods d used here, 
along with other relevant experimental parameters, are 
Hsted in Table I [H]. 
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d (/im) 


1.80(2) 


3.5(1) 


6.5(1) 


9.3(1) 


Uo (En) 


33(1) 


26(2) 


32(3) 


29(3) 




12(2) 


14(5) 


4(1) 


5(3) 


Dtf (mhi) 


55(6) 


57(6) 


45(4) 


46(5) 


DTp/d 


31(3) 


16(2) 


7(1) 


5.0(5) 



TABLE I: Experimental parameters for the four lattice peri- 
ods investigated. The Thomas-Fermi diameters Dtf along z 
were calculated using the measured atom numbers and known 
scattering length and trap frequencies. Lattice depths were 
obtained by measuring the maximum kinetic energy of an 
atom during its evolution in a lattice. The number of occu- 
pied lattice sites is approximately DTp/d. 



The lattice beams, linearly polarized perpendicular to 
the plane defined by the crossed beams, derive from a 
Ti:Sapphire laser operating at A = 810 nm (detuned be- 
low both 5S 5P transitions at 795 nm and 780 nm). 
We constructed an "accordion" lattice allowing us to con- 
tinuously vary the period of the diffracting potential [Tf] . 
Rotation of a galvanometer-controlled mirror causes the 
relative angle of the two beams to change while main- 
taining their intersection at the BEC. The radius 
of each beam is « 200 iim. The lattice is turned on 
abruptly (^ 500 ns) to its nominal depth of 30 E^, held 
constant for a variable time Tpuisc and then turned off 
abruptly. This constitutes a "pulse" of the lattice poten- 
tial. Immediately after the pulse, we release the atoms 
by turning off the magnetic trap in « 250 /US. After 
the atom cloud expands for 20.2 ms, we record the spa- 
tial distribution of the atoms using resonant absorption 
imaging |18j . Each image approximates the momentum 
distribution at the time of release. Figure 1 shows a con- 
catenated series of such images as a function of Tpuise at 
four different lattice periods. Together, the images re- 
veal the evolving momentum distribution for each lattice 
period, which collapse and revive with the characteristic 
features predicted in Ref. [S]. 



Results and interpretation 

Figure la depicts the measured momentum distribu- 
tion as a function of pulse duration for a BEC diffracted 
by a 1.8 fim lattice. For small Tpuiso, the position of 
the apparent edge of the momentum distribution grows 
linearly with Tpuisc- As expected, the distribution is 
composed of diffraction orders separated by 2hKi^. At 
early times, our data are consistent with the diffraction 
predicted using the Raman-Nath approximation. The 
Raman-Nath approximation can be viewed in a number 
of ways. For example, it neglects the kinetic energy term 
in the single-particle Hamiltonian during application of 
the pulse. This is equivalent to assuming that the only ef- 
fect of the pulse is to impose a spatially periodic phase on 
the atomic wave function, with no effect on its amplitude 
profile. This implies that the atoms move by a distance 
small compared to the lattice period d during the pulse. 
The Raman-Nath approximation is valid when 

Tpuisc ^ ^RN = „ — (2) 

where Ej^ = h^Kf^/2M is the lattice recoil energy and 
Tho = 271 /uho is the harmonic oscillator period. (Note 
that the lattice recoil in general differs from the photon 
recoil, but becomes equal to it when the lattice beams are 
counterpropagating.) In this approximation, the frac- 
tional atomic population of the n**^ diffraction order is 
J^(C/o^puise/2?i) where the J„ are Bessel functions of the 
first kind. 

Figure la shows that as the pulse duration increases 
beyond ^rn, the apparent edge of the momentum dis- 
tribution is bounded by a maximum momentum Tifcmax- 
We use this observed value of Tifcmax to determine the lat- 
tice depth Uq = h'^k'^^^/2M. The numerical calculations 
verify the accuracy of this identification to within 3 E^i- 
(The validity of the identification relies on the lattice be- 
ing deep enough to support many bound states and the 
rise and fall times of the pulse being short compared to 
Tho.) Once the edge of the distribution reaches fcmax, it 
gradually fades and a new outgoing edge appears near 
k — 0, shortly after Tho/2. Diffraction orders re-emerge 
as the new edge moves outward to higher momentum and 
the process approximately repeats. At each collapse, a 
large fraction of the population returns to the lowest or- 
ders. Collapses repeat at times which are approximate 
multiples of Tho/2 . 

The collapses and revivals can be understood from a 
classical model where atoms initially at rest in the sinu- 
soidal potential oscillate; those starting near the bottom 
of the potential oscillate at approximately Who- As we 
will see below, the suddenness of each collapse results 
from the anharmonicity of the potential and can also be 
understood in our classical model. Because of the anhar- 
monicity, no collapse is total, and there remains a sizable 
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FIG. 1: Concatenated absorption images of diffraction patterns, showing the evolving momentum distribution at four lattice 
periods d. The momentum is scaled by fcmax = 2maxM/tTOF?i where jZmax is the measured maximum amplitude of the diffraction 
pattern after time-of-fhght. The calculated initial Thomas-Fermi radius in the lattice direction, z, ranges from 26 /.im to 29 /xm, 
and in the transverse directions, from 6 /im to 10 /im. The simulations were convolved with a Gaussian to model the effects of 
finite cloud size, mean-field-driven expansion and finite optical resolution. taN is given by Eq. 2 and is the time at which the 
Raman-Nath approximation is expected to substantially fail. 
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FIG. 2: Splitting of diffraction orders (1.8 fim lattice period) 
for pulse durations beyond the Raman-Nath regime. The ef- 
fect is more pronounced the longer the pulse duration and two 
examples are indicated with white boxes. 

occupation of the higher momentum orders at the col- 
lapse point. This collapse and revival always occurs at 
times greater than and illustrates a complete break- 
down of the Raman-Nath approximation. 

Figures Ib-d depict the momentum evolution for longer 
lattice periods. While each evolution is similar, there are 
differences. For d > 1.8/im, the diffraction orders over- 
lap; indeed, mean-field-driven expansion of the individ- 
ual orders can lead to such overlap for all times-of- flight. 
The apparent discreteness in the momentum distribution 
in Fig lb, does not represent individual momentum or- 
ders but rather a modulated momentum-envelope over 
several unresolved orders. As the lattice period d in- 
creases and the spacing between the momentum orders 
decreases, the momentum distribution appears as a con- 
tinuos curve. Also, the point of collapse occurs increas- 
ingly close to Tiio/2 as d increases. 

The quantum simulations shown in the right-hand col- 
umn of Fig. 1 reproduce these features of our experi- 
ment. They depict the in- situ momentum distribution 
after Tpuisc- We numerically solve the time-dependent 
Gross-Pit aevskii (GP) equation during Tpuiso including 
the mean-field interactions [T5]. We assume that the so- 
lution factors into a time-independent radial wavefunc- 
tion and a time-dependent axial wavefunction. We treat 
the axial component using a ID GP equation with ef- 
fective interaction strength guy — Ag^n) / ii: R^Ry where 
<73D = AiTTi^as/M is the 3D interaction strength pU] . 
Here i?x and Ry are the Thomas-Fermi radii in the di- 
rections perpendicular to the lattice, and is the s-wave 
scattering length for our ^^Rb atoms. Figure 1 shows the 
good agreement between these simulations and our ex- 
perimental data. The simulation stops at the end of each 
pulse and we plot the momentum distribution (before the 
time-of-flight). 

A detail of the experimental results, absent in the sim- 



ulations, is the visible splitting of some diffraction orders 
for pulse durations beyond the Raman-Nath regime (see 
Fig. 2). The effect becomes more pronounced the longer 
the pulse duration. The effect implies spatial structure 
larger than the lattice spacing and it becomes more pro- 
nounced the longer the pulse duration. One possible ex- 
planation is that the approximate periodic translational 
symmetry of the lattice is gradually compromised as the 
inhomogeneous mean-field interaction increases over time 
due to the atoms accelerating radially toward the cloud 
center because of the dipole force arising from the laser 
beam profile. If this were the mechanism, it would ex- 
plain why we do not see the order splittings in the simu- 
lation since radial dynamics are neglected. 

Certain salient features reproduced by the numerical 
simulations can also be understood on the basis of simple 
arguments, both quantum and classical. In the following 
two sections, we first give a single-particle quantum me- 
chanical argument which explains the collapse and revival 
periods and their deviations from Tho/2. We continue 
with a classical explanation for the sudden collapse of the 
higher momentum orders and their subsequent revival. 

Quantum Mechanical Analysis 

To understand the atomic evolution after sudden ap- 
plication of an optical lattice, it is useful to decompose 
the initial BEG state into the relevant eigenstates of the 
lattice potential with energies En- We assume an ini- 
tial state with momentum eigenvalue zero, equivalent to 
an infinite-extent BEG, projected onto the Bloch states 
(with zero quasimomentum) of ID sinusoidal lattices hav- 
ing various periods. Only even parity Bloch states are 
occupied since the initial wavefunction is symmetric. We 
calculate the projection onto all significantly occupied 
states and find that the bound states contain the vast 
majority of the population. Figure 3 shows the popula- 
tions of all Bloch states up to the first unbound state, 
which include more than 99.5% of the population. 

The number of significantly occupied eigenstates grows 
with increasing lattice period. The increase of occu- 
pied bound states corresponds to the system becoming 
increasingly classical. In the deep-lattice limit, the num- 
ber of bound states in a sinusoidal potential scales as 
{Uq/ Ei^Y^"^ . Since the lattice recoil energy oc 
the number of bound states in a lattice grows linearly 
with the lattice period at fixed total depth. For example, 
a lattice with a depth of 30 En formed by two counter- 
propagating beams has a period d = A/2, a depth of 
30 i?L, and 4 bound states. However, a 30 i?R lattice 
formed by two beams intersecting at 87 mrad « 5° has a 
period d — 9.3 ^m, a depth of 15.8 x 10^ i?L, and « 80 
bound states. 

The phase of the v}^ Bloch state evolves independently 
as LOnt where aj„ = En/Ti. In a harmonic potential, the 
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FIG. 3: Calculated projections of a uniform wavefunction 
onto even-parity, zero-quasimomentum Bloch states of a 
30 -Er, lattice with periods of 1.8 /^m, 3.5 /im, 6.5 ^m, and 9.3 
fira. Each point corresponds to one such Bloch state. The hor- 
izontal axis gives the energy separation to that state's even- 
parity, higher energy neighbor, normalized to the separation 
of the lowest energy pair, and the vertical axis gives the frac- 
tional occupation of the state. The difference between con- 
secutive energy levels first decreases with increasing n, then 
increases as states become unbound. This explains the sharp 
kinks at the last data point (the first unbound state). 



frequency differences would be multiples of the harmonic 
frequency and all even eigenstates would rephase with 
period Tho/2. The anharmonic potential of a single lat- 
tice site leads to a non-uniform spacing of the energy 
levels, so the system never perfectly rephases. Neverthe- 
less, the wavefunction approximately rephases such that 
a large fraction of the population is in diffraction orders 
close to zero momentum at times close to Tho/2 . (See 
Fig. 4 where a wide spread in position corresponds to the 
collapse of the momentum distribution as seen in Fig. 1.) 



Classical Analysis 

Many aspects of this quantum mechanical system can 
be understood classically, in some cases quantitatively. 
We study the corresponding classical system by calcu- 
lating the trajectories of an ensemble of particles ini- 
tially at rest and distributed uniformly in a sinusoidal 
potential (see Fig. |5|. Each trajectory reaches a turning 
point where the velocity returns to zero. For small ampli- 
tude oscillations in the sinusoidal potential, the motion 
is nearly harmonic and atoms reach these turning points 
approximately concurrently. Larger amplitude oscilla- 
tions are increasingly anharmonic, leading to increasing 
times to reach the turning points. Figure |5] displays sev- 
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FIG. 4: Calculated evolution of the probability density 
\^{z,t)\^ within a single site of a lattice (potential shown 
in red) with period d — 1.8 /im and depth 30 Eb.- Posi- 
tions z/d = —0.5 and 0.5 correspond to the edges of the well. 
When Tpuise/Tko ~ 0.25 there is a large increase in density at 
the center of the well, at which time a large number of high 
momentum states are occupied. 



eral momentum-trajectories of initially motionless classi- 
cal particles evolving in a sinusoidal potential [21 . 

Insight can be gained from this classical picture. For 
example, the suddenness of the collapse can be under- 
stood [7^ by referring to a single-well phase space por- 
trait plotting position vs. momentum as shown in Fig. [6] 
Starting with a uniform distribution at rest {p^ — 0), the 
distribution rotates clockwise about the origin. At the 
origin, the period of rotation is 27r/who- Points farther 
from the origin rotate more slowly about the origin. The 
rotational period diverges for atoms nearest the top of the 
sinusoidal potential, z = zLd/2. The classical analog of 
our measured momentum distribution is the projection of 
this evolving phase space distribution onto the momen- 
tum axis. The horizontal tangents of the phase space 
distribution (indicated in the figure) project to singular- 
ities in the momentum distribution. The locus of these 
tangent points is referred to as a caustic. As the distribu- 
tion approaches the first turning point, the two horizontal 
tangents are near the maximum momentum; there are no 
singularities near = 0. When the central part of the 
distribution crosses the turning point (at rho/2), another 
pair of caustics suddenly emerges from the origin. This 
explains the asymmetry in the momentum evolution; the 
caustics only appear after rho/2. (The time evolution is 
symmetric only for a perfectly harmonic potential.) 

In a quantum or wave-optics system, diffraction soft- 
ens the divergence of the underlying classical caustics. 
The sequence in Fig. la-d shows, for the first time in a 
matter wave system, a progression from diffractive caus- 
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FIG. 5: Classical motion of particles in a sinusoidal poten- 
tial. The curves correspond to the trajectories of particles 
starting at different points in the potential, each with zero 
initial momentum, = 0. The smallest amplitude oscilla- 
tions correspond to the smallest period, Tho. The oscillation 
period diverges for particles starting increasingly near the top 
of the sinusoidal potential. The blue curve indicates the first 
caustic, a classical analog to the first apparent high density 
feature in our quantum system. 
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FIG. 6: Classical phase space evolution of atoms in a single 
well of a sinusoidal potential. Each black curve depicts the 
phase space distribution at a specific time. The lips of a 
single well are located at positions z = ±d/2. The phase 
space evolution is characterized by fixed points at the lips 
and rotation approaching the harmonic frequency near z = Q. 



tics [5ni EH IS2 ES] toward classical caustics: the difFrac- 
tive structure is manifest for the shortest period lattice 
and nearly invisible in the longest. 



Conclusion 

We measured the collapse and revival of the diffraction 
pattern of a BEC exposed to a pulsed ID optical lattice. 
Our results are found to be in good agreement with the 
predictions of the time-dependent Gross-Pitaevskii equa- 
tion. In addition, we employed a classical model that 
captures essential features and adds physical insight to 
the evolving diffraction pattern. For long lattice periods, 
bound states proliferate and we observed classical behav- 
ior in the long-time evolution of the momentum distribu- 
tion. We captured, for the first time, the emergence of 
both diffractive and classical caustics in ultra-cold atom 
systems interacting with optical lattices. 
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